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随 着 航天 活动 的 日 益 增 多 ， 高 轨道 飞行 器 搭载 GNSS 
点 。 以 地 球 同 步 轨 道 卫 星 TTS-2 为 例 ， 分 析 了 TTS-2 


7. 中 国 科学 院 大 学 ， 北 京 100049) 


灾 收 机 为 其 提供 导航 服务 成 为 新 的 
星 载 GPS 单 频 C1/L1 测量 数据 及 


Fo GEO 星 载 接收 机 可 接收 到 地 球 未 遮掩 的 GPS 主 久 信号 及 旁 办 信号 。C1/L1 噪声 分 


KI 


范围 较 大 ， 且 随 高 度 角 存 在 一 定 的 相关 性 。 采 用 先 验 动力 学 参数 误差 约束 的 轨道 解 算 方法 ， 利 


的 动力 学 参数 约束 策略 进行 试验 。 试 验 结果 表明 : 星 载 GPS C1/L1 数据 的 GEO 卫星 轨 
一 致 性 精度 可 达到 米 级 甚至 亚 米 级 。 在 约束 运动 状态 和 太阳 光 压 参数 的 情况 下 ，36 h 和 


(RE PUN PL EB HAE RMS 误差 均 小 于 0.65 m IX 


约束 太阳 光 压 参数 时 ，72 h IKE 


i ER- RMS REN 10m, 36h 弧 长 定 轨 的 轨道 重 车 三 维 RMS 误差 小 于 1.5 m 
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从 动力 学 参数 误差 约束 的 轨道 解 算 方 法 ， 合 理 设置 动力 学 参数 误差 ， 可 获得 更 优 的 轨道 


ial: 轨道 测定 GEO; 先 验 误差 约束 ; 漏 信号 
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随 着 全 球 卫 星 导航 系统 (global navigation satellite system, GNSS) 的 快速 发 展 ，GNSS 


逐步 成 为 “天 空地 ”用 户 获取 位 置 服务 的 重要 手段 。 大 多 数 航天 器 搭载 GNSS 接收 设备 为 
其 提 


其 提供 实时 位 置 跟踪 ， 或 提供 事后 高 精度 轨道 文 撑 其 工作 


[2,3] 
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中 于 星 载 GPS 或 Multi-GNSS 的 低 轨道 飞行 器 精密 定 轨 ， 
卫星 ”"、 风 云 三 号 卫星 ””、 天 平一 号 卫星 "等 。 近 年 来 ， 
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任务 。 然 而 ， 国 内 研究 人 员 多 数 集 
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等 高 轨道 卫星 任务 以 及 深 空 探测 任务 日 益 增 多 ，GNSS 星座 之 上 的 高 轨道 飞行 器 采用 GNSS 
进行 导航 位 置 服务 成 为 研究 热点 。 高 轨道 飞行 器 采用 GNSS 提供 自主 导航 位 置 服务 ， 不 仅 
可 以 有 效 节约 地 面 测控 资源 ， 而 且 可 以 提升 轨道 测定 精度 及 飞行 器 的 业务 能 力 。 

自 20 世纪 末 ， 欧 美 率先 开展 GPS 接收 机 在 高 轨道 飞行 器 上 的 在 轨 信 号 测试 ”， 于 21 
世纪 初 进一步 实现 了 在 轨 实 时 导航 及 定 轨 ””。 国 内 开展 高 轨道 GNSS 在 轨 测 试 研究 较 晚 。 
2014 年 ， 再 入 返回 的 嫦娥 五 号 试验 器 (Chang’E-5 Tester, CE-5T) 搭载 GNSS 接收 机 收集 了 
距离 地 心 (1~6)X104 km 的 GNSS 信号 ， 实 现 了 百 米 级 的 导航 定 轨 ”“。2017 年 ， 通 信 试 
验 技术 二 号 卫星 (No.2 Telecommunication Technology Test Satellite, TTS-2) 搭载 GNSS 接 
收 机 ， 仅 采用 伪 距 测量 和 动力 学 滤波 定 轨 ， 实 现 了 位 置 优 于 30 m 的 实时 自主 导航 “。 
内 研究 人 员 对 TTS-2 LEWER GPS 数据 定 轨 处 理 ， 采 用 重 登 弧 段 进行 比较 分 析 ， 三 维 均 
方 根 RMS (root mean square) 误差 接近 2 m. TIEA T EAA T gh GNSS 特点 
和 轨道 确定 技术 及 其 发 展 历程 ， 阐 述 了 高 轨道 GNSS 的 应 用 发 展 方向 。 

对 于 地 球 同 步 轨道 (geostationary earth orbit, GEO) 卫星 的 星 载 GNSS 接收 机 ， 可 见 
GPS 卫星 的 几何 分 布 较 差 ， 接 收 信号 功率 较 低 ， 所 处 的 空间 环境 严峻 ” 。 基 于 此 ， 本 文采 
= | TTS-2 卫星 的 星 载 GPS 卫星 数据 进行 精密 定 轨 研 究 。 首 先 介 绍 了 GEO ÆR GPS 接收 机 
的 实际 观测 条 件 ， 并 分 析 了 数据 噪声 及 其 精度 衰减 因子 (dilution of precision, DOP)。 人 然后 ， 
介绍 了 GEO 星 载 C1/L1 单 频 数据 的 测量 方程 及 其 误差 修正 。 最 后 ， 针 对 GEO 星 载 GPS 
数据 特点 ， 采 用 先 验 动力 学 参数 误差 约束 的 轨道 解 算 方法 设计 试验 ， 评 估 动 力学 参数 不 同 约 
束 策略 下 的 定 轨 重 对 一 臻 性。 本文 的 研究 内 容 对 高 轨道 飞行 器 基于 星 载 GNSS 的 高 精度 轨 
道 确定 有 重要 的 参考 价值 。 


2 GEO ÆR GPS 观测 数据 


2.1 GPS 漏 信号 的 测量 数据 

TTS-2 是 一 颗 GEO 卫星 ， 搭 载 的 GPS 接收 机 为 其 提供 实时 导航 服务 。GPS 卫星 的 
导航 天 线 指向 地 球 ， 导 航 信号 的 波束 角 约 为 42.6%”"。T 了 TS-2 运行 高 度 超过 GPS 星座 的 运 
行 高 度 ， 因 此 接收 的 GPS 信号 主要 有 未 被 地 球 遮 掩 到 的 主办 信号 和 旁 闪 信号 ， 即 漏 信号 。 
TTS-2 对 GPS 卫星 的 可 见 性 如 图 1 所 示 。 

本 文 收集 到 TTS-2 ER GPS 观测 数据 的 起 止 时 间 为 : 2018 年 的 年 积 日 (day of year, 
DOY) 100 16:00 一 108 16:00. TTS-2 星 载 GPS 观测 数据 为 C1/L1 单 频 数据 ， 数 据 采 样 间 
KAY 2 so 

TTS-2 对 GPS 星座 的 可 见 性 天 底 图 如 图 2 所 示 ( 星 固 坐标 系 ，z 指向 卫星 的 飞行 方向 ， 
z 指向 地 心 ，y 指向 zz 面 法 向 ， 构 成 右手 坐标 系 )。 接 收 机 接收 GPS 信号 的 高 度 角 最 高 约 
为 -51?。 由 于 能 接收 到 GPS 卫星 的 最 低 高 度 角 约 为 -83"*， 所 以 图 2 中 心 部 分 星空 白 。 图 
2 中 ， 方 位 角 为 245"、 高 度 角 为 -64° 和 —73° 的 两 块 区域 附 近 旦 空白， 这 是 因为 GPS 信和 号 
的 主办 与 旁 注 之 间 不 存在 导航 信号 。 
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y(Az=90) 


-x(Az=180) 


1 TTS-2 可 见 的 GPS 信号 2 ”TTS-2 对 GPS 星座 的 可 见 性 


李冰 等 人 ”指出 ，TTS-2 接收 机 天 线 在 较 大 的 波束 内 具有 较 高 的 增益 。 天 线 在 装 星 前 
对 不 同 角度 的 增益 实测 值 为 : GPS 的 L1 信号 在 —30° 的 增益 为 7.6 dB; Æ 0° 的 增益 为 
10.8 dB， 在 30° 的 增益 为 7.3 dB。 
2.2 GPS 漏 信 号 的 测量 噪声 

测量 数据 的 精度 对 精密 测定 轨 的 结果 具有 重要 影响 。 星 载 GPS 观测 数据 为 C1/L1 单 频 
数据 ， 本 文采 用 历 元 间 三 次 差 法 计算 伪 距 、 相 位 测量 值 的 噪声 ”"， 以 此 得 到 观测 数据 的 

TTS-2 平均 每 天 观测 值 的 个 数 与 高 度 角 、 噪 声 的 分 布 关 系 如 图 3 所 示 。TTS-2 可 见 
GPS 信和 号 的 高 度 角 范围 约 为 -83。~ 一 51?。 在 高 度 角 为 -81.5° ~ -79.5"， 伪 距 /相位 测量 
的 噪声 分 布 范围 明显 降低 。 整 体 上 ， 伪 距 测 量 噪声 多 分 布 于 土 40 m 之 内 ， 相 位 测量 噪声 多 
分 布 于 土 9 cm。 对 于 伪 距 测量 ， 噪 声 分 布 较 均 匀 ， 在 土 10 m 的 噪声 范围 内 存在 集中 ， 在 高 
度 角 -80"*、 噪 声 0 m 附近 存在 明显 的 集中 点 ， 观 测 数据 的 个 数 达 到 8X103 左右 。 在 高 度 角 
—70° 和 高 度 角 -62。 附近 存在 轻微 的 集中 。 对 于 相位 测量 ， 在 土 3 cm 的 噪声 范围 内 集中 。 
在 高 度 角 —78°, —70° 和 —62° 附近 存在 明显 的 集中 ， 观 测 数据 的 个 数 大 于 1.8X104。 经 统 
计 ，DOY 100 16:00 一 108 16:00 的 伪 距 测量 值 噪声 RMS 误差 为 7.3 m， 相 位 测量 值 噪声 
RMS 误差 为 1.3 cm. 

以 G29 卫星 为 例 ， 我 们 研究 了 伪 距 /相位 测量 值 噪声 随 高 度 角 的 变化 关系 ， 结 果 如 图 4 
所 示 。DOY 101 5:00— 6:00 的 伪 距 测量 值 噪声 变化 区 间 为 (—20, 20) m， 相 位 测量 值 噪声 
变化 区 间 为 (—4, 4) cm， 高 度 角 自 —67° 降低 至 -73"。 自 6:30 左右 开始 ，G29 卫星 再 次 被 
TTS-2 卫星 可 见 ， 高 度 角 从 -75。 逐渐 下 降 至 -80*， 随 后 逐渐 升 高 至 -66*。 在 高 度 角 升 高 
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图 3 


高 度 角 /(*) 


TTS-2 平均 每 天 观测 值 个 数 与 高 度 角 、 噪 声 的 分 布 关 系 


2018 DOY 101 时 刻 /h 


图 4 G29 PẸ C1/L1 噪声 随 高 度 角 的 变化 
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的 过 程 中 ， 其 变化 速率 越 来 越 绥 。 在 此 期 间 ，C1/L1l 测量 值 噪声 的 分 布 范围 随 高 度 角 变化 
而 变化 。 当 高 度 角 小 于 -75。 IN, C1/L1 测量 值 的 分 布 范 围 分 别 是 土 3 m 和 土 2.2 cm。 随 着 
高 度 角 的 升 高 ，C1/L1 测量 值 噪声 的 分 布 范 围 逐 步 扩 大 。 在 G29 卫星 不 可 见 前 的 一 小 段 时 
间 ，C1/Z1 测量 值 噪声 快速 变 大 。 

2.3 TTS-2 对 GPS 的 可 见 星 和 DOP 

GPS 测量 中 ， 接 收 机 单 历 元 的 可 见 卫 星 数 及 其 分 布 直接 影响 位 置 服务 精度 。 位 置 精度 
衰减 因子 是 评判 GPS 导航 定位 精度 水 平 的 参考 指标 。 对 于 GEO ÆR GPS 测量 ， 水 平 精 度 
衰减 因子 (horizontal dilution of precision, HDOP) 和 垂直 精度 衰减 因子 (vertical dilution of 
precision, VDOP) 分 别 表示 GEO 卫星 轨道 迹 向 与 法 向 平面 的 精度 衰减 因子 和 径 向 精度 衰减 
因子 。 较 低 的 精度 衰减 因子 说 明 可 获得 较 高 精度 的 几何 位 置 解 。 

TTS-2 星 载 接收 机 单 历 元 的 GPS 可 见 星 数 、HDOP 和 VDOP 如 图 5 所 示 。 接 收 机 仅 
设计 了 8 个 信和 号 通道 ， 最 高 可 接收 8 MI GPS 卫星 。GEO 星 载 接收 机 在 大 部 分 历 元 对 GPS 
可 见 星 数 为 6 ~ 8 颗 ， 少 部 分 历 元 对 GPS 可 见 星 数 为 4 颗 或 5 颗 。 大 部 分 历 元 的 HDOP 值 
集中 于 2 ~ 10，VDOP 值 集 中 于 6 ~ 30。 


可 见 卫星 数 
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5 TTS-2 星 载 接收 机 的 GPS 可 见 星 数 、HDOP 和 VDOP 


Chinay i 
CNINAA Ivi 


122 天 文学 进展 40 卷 


3 GPS 测量 方程 与 误差 修正 


31 ”测量 方程 
GEO ÆR GPS 接收 机 的 C1/L1 测量 方程 ， 不 考虑 硬件 延迟 ， 可 以 分 别 表示 为 : 


LI = p +c: (dt, — dt) — Ay,» NB 


一 7o + Trop? + eri 


er = p> +c: (dt, — dt®) + Ion’ + Trop’ +e, , 


AH, 表示 GPS 卫星 s 天 线 中 心 到 接收 机 r 的 理论 距离 ，dt. 和 de® 分 别 表示 接收 机 r+ 和 
GPS 卫星 s 的 钟 差 ，c 表示 真空 中 的 光速 ，Jons 和 Tro 分 别 表示 信和 号 传播 过 程 中 的 电离 层 
延迟 和 对 流 层 延迟 ，s2; 和 6, 分 别 表 示 GPS 卫星 s 伪 距 C1 和 相位 L1 测量 值 的 噪声 ， 和 yp 
ALN 分 别 表 示 频 率 fa 真空 中 传播 的 波长 和 未 知 的 模糊 度 ，C1s 和 L 分 别 表示 接收 机 
测定 的 GPS 卫星 s 的 伪 距 C1 和 相位 L1 测量 值 。 
3.2 ”硬件 延迟 修正 

采用 国际 GNSS 分 析 中 心 (IGS) 提供 的 事后 精密 轨道 和 钟 差 产品 ， 对 TTS-2 星 载 GPS 
单 频 C1/L1 测量 数据 进行 处 理 。IGS 产品 采用 双 频 P1/P2 伪 距 和 L1/L2 相位 的 无 电离 层 
组 合 观 测量 解 算 ， 因 此 在 处 理 单 频 C1/L1 测量 数据 时 需要 考虑 不 同 频 率 信 号 的 硬件 延迟 。 
单 频 相位 测量 与 双 频 相位 测量 的 硬件 延迟 时 间 差 可 被 模糊 度 和 接收 机 钟 差 吸收 ， 故 本 文 不 再 
考虑 相位 测量 的 硬件 延迟 时 间 差 。 采 用 欧洲 定 轨 中 心 (the Center for Orbit Determination in 
Europe, CODE) 的 差分 码 偏差 (different code biase, DCB) 产品 修正 C1 伪 距 与 双 频 P1/P2 
伪 距 的 硬件 延迟 时 间 差 ， 其 修正 方法 ”表示 为 : 
DCBP1p2 

fr- 

AF, DCB, p. 和 DCBS,¢, 表示 GPS FÆ s E P1 分 别 与 伪 距 P2 和 伪 距 C1 的 差分 
码 偏差 ，Cl9cp 表示 GPS TÆ s 经 过 DCB 修正 后 的 伪 距 。Cl5cp 可 以 直接 采用 IGS 精密 
产品 进行 处 理 。 
3.3 ”电离 层 与 对 流 层 延迟 修正 
对 流 层 分 布 于 地 表 以 上 20 km 以 下 ， 电 离 层 分 布 于 地 球 上 空 50 ~ 1000 km。 如 图 1 所 
I> GEO 星 载 接收 机 在 高 度 角 较 低 时 接收 到 的 GPS 信和 号， 贴近 地球 表面 经 过 电离 层 和 对 流 
层 ; 当 高 度 角 较 高 时 ，GPS 信和 号 不 经 过 电离 层 和 对 流 层 而 直接 被 接收 。 设 GEO 卫星 的 轨道 
半径 为 Rano GEO 卫星 接收 到 的 GPS 信号 传播 路 径 距 地 心 的 最 近 距 离 为 h， 则 接收 信号 
的 高 度 角 ele 表示 为 : 


Clbceg = Cl’ +c: ( fi + DCB) ， (2) 


h 
ele = arctan — 90° . (3) 
Raro 


设 地 球 半径 为 Rear “4 (hh 一 Rearth) 大 于 电离 层 高 度 时 ，GEO 接收 到 的 导航 信号 不 受 


电离 层 和 对 流 层 的 影响 。 本 文 参与 轨道 解 算 的 GPS 测量 数据 的 截止 高 度 角 设 为 —78°, AR 
经 过 电离 层 和 对 流 层 的 导航 信号 测量 数据 。 硬 件 延 迟 修 正 后 ， 忽 略 电离 层 和 对 流 层 延迟 ， 式 
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(1) 可 以 写作 : 
| 


(4) 
L15 = o +c- (dt, — OF) Ape NG teha « 
A 先 验 动力 学 参数 误差 约束 的 轨道 解 算数 学 模型 
卫星 在 时 刻 t 的 运动 状态 可 以 表示 为 : 
= G (Tto Vio» f) ) (5) 
Vr 


RP, a, v RE EBE t NAIM MAME, £o, v, 表示 卫星 在 to 时 刻 的 位 置 和 速度 ，f 
表示 待 估 摄 动力 参数 ，G@G 表示 卫星 在 to 时 刻 运动 状态 向 t 时 刻 运动 状态 的 映射 函数 。 式 (4) 
中 的 理论 距离 pos 可 表示 为 p= H(a,, 0°). H a GPS TEME x Al GEO TEME a, 
到 理论 距离 的 映射 函数 。 对 于 时 刻 t WAE LOs (t) 和 相位 LL (t) 可 以 表示 为 : 


a (t) = HIG(zi, » Utos f), x, dt,] + ELC1(t) ; (6) 


LAY (t) = H|G (£to, Vto; T) x, dtr, 和] 十 ELE) ， 


RP, ero, erne 分 别 表示 观测 数据 伪 距 LO ben (t) 和 相位 LE (t) 测量 的 噪声 。 一 般 解 
算 轨道 时 ， 首 先 给 定 eu, ws 和 上 的 近似 值 ， 代 入 式 (6) 中 进行 迭代 求解 。 

GEO 卫星 的 运动 状态 受到 动力 学 约束 ， 动 力学 参数 的 变化 在 短 时 间 内 非常 小 。 同 时 ， 
动力 学 定 轨 的 法 方程 是 病态 方程 ， 容 易 受到 观测 量 的 粗 差 影响 ， 使 方程 的 解 无 法 达到 最 优 。 
然而 ，GEO BAX GPS 的 测量 数据 噪声 较 大 ， 部 分 数据 包含 粗 差 ， 且 每 颗 星 噪声 的 分 布 范围 
不 稳定 ， 容 易 对 定 轨 解 算 的 病态 法 方程 的 解 偏离 最 优 值 。 因 此 本 文采 用 先 验 动力 学 参数 误差 
约束 的 轨道 解 算 方法 ， 其 主要 思路 如 下 。 

将 轨道 解 算 的 动力 学 参数 初 值 zw， wx 和 了 视 为 虚拟 观测 量 ， 并 根据 其 先 验 误差 赋予 相 
应 的 权重 参与 轨道 解 算 ， 有 利于 定 轨 处 理 中 剔除 粗 差 和 噪声 较 大 的 观测 量 ， 使 法 方程 解 得 最 
优 解 。 对 于 式 (6) 组 成 的 观测 方程 组 ， 增 加 初始 动力 学 参数 zu, v1, 和 /的 虚拟 观测 量 : 


Tto Lz,, Var, 
Uto = Ly,, + Vor, ) (7) 
f Ly Vy 


AP Le, Lons Le M Ve 


( 即 初 值 ) 及 改正 值 。 
模型 的 观测 向 量 为 : 


to? Vono Vp 分 别 表示 初始 动力 学 参数 Lio, v1。 和 f 的 虚拟 观测 量 


T 
[z Ly, Lj, LODosg (ti), ,了 Clpce (ta), LEF (t,),--- | 


Tto? 
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观测 向 量 的 权 阵 为 : 


diag 2 ) 7 1-3 3 SS oe | 4 st 


2 2 X ad 
ILa ILo OL; Probe (t) PTO beB (tn) OLr1s (t1) OLEIS (tn) 


5 定 轨 策 略 与 结果 分 析 


5.1 ” 定 轨 策略 与 试验 设计 

本 文采 用 先 验 动力 学 参数 误差 约束 的 轨道 解 算 方法 ， 对 动力 学 参数 采用 不 同 的 约束 策 
略 进行 轨道 处 理 ， 评 佑 定 轨 精度 。 设 计 了 三 个 试验 ， 具 体 分 析 过 程 如 下 : 
(1) 约束 动力 学 参数 为 卫星 的 初始 位 置 、 速 度 和 5 个 光 压 参数 。 以 前 一 次 邻近 弧 段 的 定 
轨 解 算 参 数 及 其 中 误差 作为 本 次 定 轨 参数 初 值 和 误差 rc; 所 有 待 估 参 数 的 误差 分 别 采用 lc 
或 3c， 试 验 分 别 记 为 Test-1a 或 Test-1b. 

(2) 约束 动力 学 参数 为 5 个 光 压 参数 。 以 前 一 次 邻近 弧 段 的 定 轨 解 算 参数 及 其 中 误差 
作为 本 次 定 轨 参数 初 值 和 误差 o; 5 个 光 压 参数 的 误差 分 别 采用 lo 或 3r， 试 验 分 别 记 为 
Test-2a 或 Test-2b。 

(3) 约束 动力 学 参数 为 光 压 的 5 个 光 压 参数 。 将 Test-la 的 各 组 太阳 光 压 参数 及 其 中 误 
差 取 平均 值 ， 作 为 光 压 参数 初 值 及 其 误差 o; 5 个 光 压 参数 的 误差 分 别 采 用 lc 或 3c， 试 验 
分 别 记 为 Test-3a 或 Test-3b。 

试验 初始 时 ， 对 试验 Test-la 第 一 组 数据 迭代 解 算 动 力学 参数 ， 获 取 其 初始 位 置 的 中 误 
差 小 于 5 cm 的 动力 学 参数 ， 作 为 三 个 试验 的 第 一 组 定 轨 待 估 动力 学 参数 的 先 验 值 。 

本 文 定 轨 所 采用 的 测量 模型 与 动力 学 模型 见 表 1。 
5.2 ”试验 结果 与 分 析 

本 文 定 轨 的 数据 起 止 时 间 为 2018 年 的 DOY 100 16:00 一 108 16:00， 分 别 进行 72 h 弧 长 
Al 36 h 弧 长 的 定 轨 解 算 。72 h MRE, BUS OAM 16:00 至 第 4 天 的 16:00， 共 可 进行 六 
组 轨道 解 算 ;36 h MK, BN OKA 16:00 至 第 3 天 的 04:00， 共 可 进行 七 组 轨道 解 算 。 
各 试验 定 轨 后 ， 伪 距 / 相 位 测量 值 残 差 及 其 数据 有 效率 统计 见 表 2。 

HX 2 可 知 ，72 h 弧 长 定 轨 ， 其 观测 数据 有 效率 为 85% AA, DEMER RMS 值 
为 537.8 cm 左右 ， 相 位 测量 值 的 RMS 值 为 5.6 cm 左右 。 各 试验 的 数据 有 效率 、C1/L1 W 
ERA RMS 值 相差 不 大 。36 h MER, Cl 测量 值 的 RMS 残 差 约 为 523.4 cm，L1 il 
量 值 的 RMS 残 差 为 5.3 cm， 观 测 数据 有 效率 为 85%。 

T2hMKNCH, BB 48 h， 共 有 五 组 重 闪 弧 段 。 每 组 重合 弧 段 比较 时 ， 为 避免 轨道 
的 边缘 效应 ， 取 重合 弧 段 中 间 的 24 h 进行 比较 。36 h MKH, ES 12 hp， 共 有 六 组 重 
BIBL, WEE 6 h 进行 比较 。72 h IHEM 36 h 弧 长 定 轨 的 各 个 试验 重合 弧 段 比 较 的 轨道 
三 维 、 径 向 、 迹 向 和 法 向 的 RMS 误差 如 图 6 和 图 7 所 示 。 
图 6 所 示 的 72 bh 弧 长 定 轨 重 共 比 较 中 ， 运 动 状态 受到 约束 的 试验 Test-la, Test-1b 的 三 
维 RMS 误差 均 小 于 1m， 径 向 RMS 误差 小 于 0.2 m。 仅 约束 太阳 光 压 参数 的 试验 Test-2a， 


= 


[fini 
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表 1 定 轨 采用 的 测量 模型 与 动力 学 模型 


模型 参数 
定 轨 弧 长 72 h/36 h 
采样 间隔 60s 
观测 量 定 权 C1: 700 cm; L1: 6 cm 
截止 高 度 角 —78° 
GPS 轨道 和 钟 差 CODE 最 终 产 品 
GPS 天 线 相 位 中 心 采用 文件 IGS14.atx 改正 
TTS-2 卫星 GPS 接收 天 线 相 位 中 心 “采用 标 称 值 ， 不 考虑 天 线 相 位 中 心 变化 
天 线 相位 中 心 缠绕 考虑 ，L1 单 频 改正 
为 距 硬 件 延 迟 修正 
电离 层 与 对 流 层 不 考虑 
模糊 度 浮 点 解 
重力 场 模 型 EGM2008 10X10 
N 体 引 力 JPL DE421 
地 球 自转 定向 参数 IERS2010 C04 
岁差 章 动 模型 IAU2000R06 
潮汐 模型 海洋 、 大 气 、 固 体 潮 、 海 潮 
光 压 模型 ECOM5 参数 模型 
KAHR 不 考虑 
经 验 力 不 考虑 


表 2 ”CI1/Z1L 测量 值 残 差 及 数据 有 效率 统计 


试验 72 h MRA 
C13Æ/cm L13% /cm 效 数 据 数 “剔除 数据 数 ”有 效率 / 史 
Test-la 536 5.63 25 862 4 508 85 
Test-1b 533 5.48 25 719 4 651 85 
Test-2a 539 5.64 25 897 4 473 85 
Test-2b 537 5.58 25 820 4 549 85 
Test-3a 540 5.69 25 953 4 417 85 
Test-3b 539 5.56 25 817 4 553 85 
试验 36 h MR EM 
C1 残 差 /cm 了 工 1 残 差 /cm ”有 效 数据 数 ”剔除 数据 数 ”有 效率 /%% 
Test-la 535 6.00 12 999 2 246 85 
Test-1b 529 5.44 13 015 2 231 85 
Test-2a 521 5.17 12 946 2 299 85 
Test-2b 523 5.29 13 036 2 209 86 
Test-3a 509 4.58 12 607 2 632 83 


Test-3b 520 5.17 12 942 2 304 85 
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图 7 


36 h 定 轨 ，12 h 重 又 比较 的 三 维 、 径 向 、 迹 向 和 法 向 的 RMS 误差 
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Test-2b, Test-3a 的 三 维 RMS RÆK, KA RMS 误差 仍 小 于 0.3 m Test-3b 的 三 维 

RMS 误差 约 为 0.65 m, K RMS 误差 小 于 0.3 m。 明 显 的 特征 为 ， 约 束 卫 星 运动 状态 

参数 的 轨道 解 算 ， 获 得 了 更 好 的 轨道 重 亚 结果 ， 且 其 轨道 法 向 误差 RMS 明显 降低 。 

图 7 所 示 的 36 h 弧 长 定 轨 重 钱 比较 中 ， 运 动 状态 受到 约束 的 试验 Test-la 的 三 维 RMS 
误差 均 小 于 0.8m, 4 RMS 误差 小 于 0.2 m; 试验 Test-lb 的 后 三 组 重合 误差 RMS RK, 
但 整体 的 三 维 RMS 误差 小 于 1 m。 仅 约束 太阳 光 压 参数 的 试验 Test-2a, Test-2b, Test-3a 
和 Test-3b 的 三 维 RMS 误差 为 0.9 ~ 1.5 m。 与 72 h 弧 长 定 轨 相 同 的 是 ， 约 束 卫 星 运动 状 
态 参 数 的 轨道 解 算 有 更 好 的 轨道 重 又 结果， 轨道 法 向 的 误差 RMS 最 低 。 

为 了 更 具体 地 评价 TTS-2 Bak GPS 单 频 数据 的 定 轨 结 果 ， 汇 总 统计 试验 Test-1, Test-2 

和 Test-3 三 个 试验 整体 的 轨道 重 受 RMS 误差 ， 见 表 3。 在 受到 卫星 运动 状态 参数 约束 的 情 
况 下 ，72 h 和 36 bh 弧 长 定 轨 均 可 获得 三 维 RMS 误差 小 于 0.65 m 的 轨道 ，72 h 弧 长 的 径 向 
RMS 误差 小 于 0.2 m, 36 h 弧 长 定 轨 的 径 向 RMS 误差 小 于 0.3 m。 在 仅 约束 太阳 光 压 参数 
的 情况 下 ， 采 用 3c 先 验 误差 约束 的 轨道 重 老 一致 性 更 好 。 综 合 而 言 ，72 h 弧 长 测定 的 轨道 

qm 重 车 一 致 性 小 于 1 m, 36 h 弧 长 测定 的 轨道 重 亚 一 致 性 小 于 1.5 m。 与 文献 [15, 16] 的 定 轨 结 

> 果 相 比 ， 轨 道 重 又 一致 性 提升 了 40% WE. 


表 3 72 h 和 36 h 弧 长 的 轨道 重 又 RMS 误差 统计 /m 

试验 72 h 弧 长 定 轨 36 h MR EM 

=4 @ Ye Bem | 三维 径 向 切 向 er 
Test-la | 0.62 0.18 0.59 0.03 | 0.63 0.29 0.56 0.03 
Test-lb | 0.63 0.15 0.59 0.12 | 0.58 024 0.5 0.13 
Test-2a | 0.96 0.23 0.82 0.36 | 145 0.56 1.15 0.5 
Test-2b | 0.77 0.17 0.64 0.35 | 1.05 0.41 0.76 0.48 
Test-3a | 0.89 0.24 0.74 0.38 | 13 047 1.05 0.51 
Test-3b | 0.65 0.15 0.53 0.3 | 0.91 0.33 0.68 0.37 
汇总 0.75 0.19 0.65 0.26 | 0.99 0.38 0.78 0.34 


6 4 论 


本 文 以 TTS-2 星 载 GPS 观测 数据 为 例 ， 对 GEO HA GPS 单 频 C1/L1 数据 进行 了 分 
析 ， 采 用 先 验 动力 学 参数 误差 约束 的 定 轨 解 算 方法 进行 试验 并 分 析 结 果 ， 主 要 结论 如 下 。 

(1) GEO 星 载 GPS 数据 噪声 分 布 范 围 较 大 ， 且 随 高 度 角 存在 一 定 的 相关 性 。 在 高 度 角 
为 一 78°, —70° 和 一 62° 附近 ， 观 测 数据 存在 较为 明显 的 聚集 效应 。GEO 星 载 接收 机 在 大 部 
分 历 元 对 GPS 可 见 星 数 为 6 ~ 8 Hi, HDOP 值 集中 于 2 ~ 10，VDOP 值 集中 于 6 ~ 30. 

(2) 星 载 GPS 单 频 数据 的 GEO La ea ea MBEAN lm 量 级 ， 最 优 可 达 
到 亚 米 级 。 约 束 卫 星 运动 状态 参数 的 轨道 解 算 ，72 h 和 36 bh WKE HIRE T ERM 
性 优 于 0.65 m 的 轨道 。 仅 约束 太阳 光 压 参数 ，72 h 弧 长 定 轨 的 轨道 重 堆 三维 RMS REN 
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1m, 36 h MR ERNE EZ =E RMS 误差 小 于 1.5 m 

(3) 试验 证 明 ， 合 理 设置 动力 学 参数 误差 ， 采 用 先 验 动力 学 参数 误差 约束 的 定 轨 解 算 方 
法 ， 可 以 得 到 重 益 一 致 性 更 优 的 轨道。 

本 文 的 定 轨 方法 、 研 究 内 容 与 结果 对 高 轨道 飞行 器 的 星 载 GNSS 高 精度 应 用 具有 一 定 
的 借鉴 意义 和 参考 价值 。 北 斗 全 球 卫 星 导航 系统 的 建立 ， 将 为 高 轨道 飞行 器 的 星 载 GNSS 
导航 及 精密 定 轨 提 供 更 好 的 解决 方案 。 另 外 ， 高 轨道 飞行 器 的 星 载 GNSS 数据 的 噪声 较 大 ， 
其 误差 来 源 仍 竺 研究。 高 轨道 飞行 器 分 米 级 甚至 厘米 级 的 高 精度 测定 轨 的 挑战 依然 存在 ， 仍 
须 进 行 深入 研究 。 
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Abstract: With the rapid development of space activities, the high-orbit spacecraft would 
carry GNSS receivers to navigate in space, which is becoming a new research focus. This 
paper takes the geosynchronous orbit satellite TTS-2 as an example, analyzes the measure- 
ment data quality and noise characteristics of the onboard GPS C1/L1 single-frequency. The 
GPS signals received by the GEO satellite are from the side lobe signals and the uncovered 
part of the GPS main lobe signals. C1/L1 noise has a large distribution range, and there is 
an obvious relativity with the elevation angle. This paper adopts priori dynamic parameter 
error constraint orbit method, carries out orbit determination experiments with different 
strategies of priori dynamic parameter constraint. The experiment results show that the 
orbital overlap consistency of GEO satellite orbit determination with GPS C1/L1 onboard 
data is about 1 m or even sub-meter level. Under the condition of constrained satellite 
motion state and solar pressure parameters, the three-dimensional RMS errors of 36 h and 
72 h arc orbit determination are both less than 0.65 m. When only constraining the solar 
pressure parameters, the three-dimensional RM S error of the 72 h arc orbit determination is 
1.0 m, and the three-dimensional RM S error of the 36 h arc orbit determination is less than 
1.5 m. The better orbit solution would be solved when the priori dynamic parameter error 
constrained orbit solution method and the reasonable dynamic parameter error restraint 
are adopted. The research methods and experiments in this paper have certain reference 


meanings for the high-precision applications of spaceborne GNSS on high-orbit spacecraft. 
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